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Summary 

The goals of the grant fall in two main areas. The first 
area Is control system design and the goals are: 1 ) to define the 
limits of "digitized S-Plane design techniques" vs. sample rate, 

2) to show the results of a "direct digital design technique" and 

3 ) to compare the two methods. The second area is to evaluate the 
roughness of autopilot designs parametrically versus sample rate. 

Goals of the first area have been addressed by 1) an 
analy. of a 2nd order example using both design methods, 2) a 
linear analysis of the complete 737 aircraft with an autoland obtained 
using the digitized S-plane technique, 3) linear analysis of a high 
frequency 737 approximateion with the autoland from a direct digital 
design technique, 4) development of a simulation for evaluation of the 
autopilots with disturbances and nonlinearities included. 

An important aspect of the results to date is that the 
damping of the short period roots of the full 737/autoland system 
degrades by 50$ at 5 cps (decreasing to 15$ at 20 cps) with a slightly 
larger degradation of the dutch roll roots when using, in our judgement, 
the best digitized S-plane design technique. 

Goals of the second area, roughness evaluation, have been 
addressed by defining an experiment to be carried out on the Langley 
motion simulator and coordinated with analysis at Stanford. 
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EXAMPLE: Pitch SAS - classical des ign 

This system under consideration for the example is shown in Figure 

t 

1. The transfer function 8 (s) was derived for a 737 in cruise mode. 

6e(s) 

The washout for nitch rate was chosen as being typical for a stability 
augmentation system. The washout break frequency, 1/T , is chosen from 
the characteristics of the gyro bias and the expected nature of aircraft 
maneuvers producing a steady state pitch rate. For this system T=l and 
K=-.35 yields roots of -2.52 ± 2.57 j and -.9865. These will be con- 
sidered the desired root locations for this example. The goal of a digital 
mechanization here will be to duplicate these roots. 

Assuming a zero-order hold (Z.O.H.) signal reconstruction, the 
z-transform of the plant is: 


tf(z) 

6e(z) 


=.5275^ 


[e 4T (cos3.28-3.24sin3.28T)-l]z+ e 


-1 4T -2 8T 

(cos3.28T+3.24sin3.28T)-e 


2 _ -1.4T _ OQ(p -2.8T 

z - 2e cos3.28Tz + e 


where T is the sampling rate. 

This is an exact result for the plant since there is a physical ZOH 
in the D/A conversion. There are, however, several possible mechanizations 
of the feedback washout. Methods under consideration and their results 
for the washout (including a constant to produce unity gain at high 
frequency) are: 
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With the exception of Tustin' s method (which also produces this 

result for high sampling rate) these methods produce the same transform, 

one that causes the pole and zero of the washout to be mapped into the 

sT 

z-plane by the transformation z=e for all sampling rates. This is 
inherently satisfying since the poles and zeroes are mapped exactly with 
no shifting at low sampling rates. 

The z-tran?forms of the plant and the washout were used to study 

the variation in the roots of the sampled system as a function of sampling 

rate (shown in Figure 2). This variation comes from the change in porition 

of the plant zero as a function of sampling rate. The feedback raachaniza- 

tion has been chosen to produce no shift in poles or zeroes and examination 

of the denominator of 6 (z) reveals poles in the z-plane corresponding to 

&e(z) 

the s-plane poles for all sampling rates. However the plant zero does not 
remain at -1.06 (the value for the continuous system), but shifts negatively 
as the sampling rate decreases (-1.1 for T=.2sec). This zero shift occurs 
for all plants with zeroes not at the origin. The variation in complex 
roots here tends to be more highly damped and the real root moves only 
slightly toward the origin. While this would seem to be an argument for 
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unconcern with slower sampling, this Is not true for all situations 
(to be seen later in the 737 analysis). 

This variation was generated with K and t being the constant values 

resulting from continuous analysis. If they are •'llowed to vary as a 

function of sampling rate a system is produced that minimizes the 

distances between the roots and the desired roots, d=yr -r 

*■*[ i idesiredl 

Figure 3 shows d for both the constant and variable gain cases, and the 
ratio of gains for the two cases. Variation of T did not cause any 
decrease in d. Varying K yielded a maximum of 40 percent improvement 
in d over the constant gain case at 7 samples/sec. 

The nature of the variation can be studied approximately by the 

following method: A root locus is drawn for the continuous system and the 

closed loop poles noted. The variation of the plant zero with sample rate 
can be computed explicitly from the expression of its z-transform. A 
family of loci may be sketched for different plant zero locations in the 
s-plane and the nature of the change in feedback may be explored. An 
Axact analysis would require z-plane root loci; however, this method gives 
the same trends. For this example, the sketch (not to scale) would be as 
in Figure 4. 

The relative invariance of the real root can be seen as well as the 
mo.ement in the complex closed loop poles shown in Figure 1. Reducing 
the magnitude of K as indicated in Figure 3 will tend to move the poles 
along the dotted lines toward the open loop poles and consequently 
toward the desired root location. 

It can be argued that the complex roots should be matched only, 
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FIGURE 3 VARIATION IN ROOT DISTANCE AND GAIN AS A FUNCTION OF SAMPLING RATE 


variation of closed 
loop root 



leaving the real root (which Is often at high frequency or relatively 
Insensitive to changes) tre' . 

There are two variables, K and T, in the characteristic equation 
and therefore they should be sufficient to specll, two roots. This is 
invalid however. If the two complex roots are 83-R^lI^j and the real 
root Is s=-R^ , the characteristic equation Is 

s 3 + (2R +R )s a + (R 3 +I 3 +2R R )s + R (R 3 +I 2 ) =0 
18 11 18 *li 

The characteristic equation as a function of K and T is of the form 

s 3 ♦ f (K,T)s 2 + f (K,T)s + f (K.T) « 0 
X 3 3 

and It can be seen by examining the dependence of R^ and 1^ on f^ , f^ , 

and f that in general K and T cannot be varied so as to keep R and I 
3 K l l 

constant . 

In conclusion, the relatively simple procedure of changing feedback 

gains as a function of sampling rate will not retain the desired closed 
loop roots although some improvement is possible. It is also not possible 
in general to match some roots, leaving the others free. The reason that 
changes are occuring in the roots for slower sample rates is that an 
additional lag is being introduced into the system. To negate the effect 
of this lag, sophistication over and above simple gain changes is re- 
quired to restore the system to its original performance. 
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EXAMPLE : State - Space Design 


Another feedback technique Is the design of an estimator and the use 
of estimated states In a control law. The s*-te-space equivalent to the 
washout Is a bias estimator tor the steady-state component in q. Figure 
5 shows the block diagram for such a system. There are two basic mech- 
anization choices for an estimator: full state and reduced state. In 

this application it is natural to choose a reduced state estimator. 


The state-space representation of V 2 system with states q and or is: 
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Computing the matrices 9 = e and f = jT e Gdt the discrete system is 


x = 


n+x 


9 x + Tu 
n n 


with the bias equation 

b = b 
n+x n 

Ignoring 6 e for the moment, a reduced state estimator can be built 
resulting in 


(1-L )<z-l) 
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q(z) 

y(z) 


' L 9 ‘ 

z -9 + 1 xa 

aa nr* 
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for gains L and L . 
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FIGURE 5 STATE SPACE MODEL 




Comparing this to the digital washout where 



the estimator transfer function contains another pole and zero. There are 
only two gains and therefore the locations of the zero and both poles are 
not independent. While an estimator adds flexibility in choosing pole 
locations, it also adds a zero whose location is determined by the poles. 
In most modern control designs the zerots between the observations and 
the estimated states are generally not taken into account because they do 
no - ,->'ar in the transfer function between the commanded and actual out- 
put. But they do appear in the transfer function between the bias and 
the output. Thus their position is more important here and due to the 
lack of flexibility in placing the z-ro, it may be difficult *-o use this 
technique for sample rate compensation. 
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Discrete Control Laws 


The continuous control laws /or the 737 autopilot in the autoland 
mode prior to flare have been digitized using the "matched z-transform"^ 1 ^ 
technique. As discussed in the example, the technique is almost identical 
to the prewarped tustin; however, it gives a better approximation when 
there are zeroes. For the transfer functions contained in the 737 auto- 
pilot, the digitized versions are: 


K 

c 

s + a 


z + e 


at 



picked to match low frequency gain 


K S 
c 

s + a 


¥ ( z— 1 ) 

i; 




picked to match hi frequency gain. 


1 z + 1 

s z - 1 


T;ie autopilot used in the digitization is shown in Figures 6, 7, 8 and 

represent our interpretation of the actual autoland autopilot. The aircraft 

linear model for the 2. 5° approach configuration has been supplied by Langley 

and was used as is. Our interpretation of the auxilary equations relating 

the autopilot output to the aircraft input is shown in Fig. 9. 

The effect of these digitizations were examined by a discrete linear 

analysis of the entire system composed of the autopilo*, aux. equ& and a/c 

model. Since there is negligible coupling between the longitudinal and 

lateral motion, the analysis was perform->d separately in the two axes. The 

z-plane roots resulting from this analysis were then transformed back to the 

lnz 

s-plane by the transformation s = — — so that they could be compared on an 
equivalent basis. 
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For the longitudinal system the continuous system (equivalent to 
Infinitely small sampling interval) has the pole locations shown in Figure 10. 
As the sampling rate decreases the only roots that move appreciably are those 
associated with the short period mode (Figure 11). A system with a short 
period plant model, 0 feedback with a washout, and an actuator lag model 
produces the family of loci of Figure 12 vs. movement of the plant zero. 

These loci are relevant because, as was discussed in the example, the plant 
-ero movement is a primary effect of discrete control. The direction of 
the movement corresponds well to the change in short period roots experienced 
by the entire longitudinal system. 

The difference between this and the example considered earlier is 
the inclusion of the actuator lag (and the fact that the open loop roots 
are different). Here however the direction of movement of the short period 
poles is approximately perpendicular to the loci, and little improvement in 
pole position will be derived from decreasing the gain based on discrete 
analysis as in the example. Using this high frequency approximation, 
modifications to the discrete washout have been examined that increase 
the short period root damping back to the desired location for the slower 
sampling rates. The modifications consist of the addition of a pole and 
zero to increase the compensation lead so as to make up for the lag due 
to the sampler. At a 5 cps sample rate with short period closed loop 
roots at 1 cps, we do not anticipate any significant noise or disturbance 
amplification due to the increased compensation lead. 
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FIGURE 7 
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FIGURE 11 SHORT PERIOD ROOT LOCATION VS. SAMPLING RATE 
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The desired lateral roots are shown in Figure 13 and the variation 

in these roots that change is shown in Figure 14. The roots that change 

5 

are thosj associated with the loop on the lateral autopilot drawing 

(Figure 7) and the roll subsidence as well as the aileron lag and the dutch 
roll. The lateral system has not yet been analyzed to determine the critical 
feedback loop where additional compensation can be introduced In a direct 
digital design to reduce root movement. These results were obtained by a 
linearization that assumed the spoiler threshold was not exceeded. In 
obtaining the final digital designs, an additional linearization where 
the threshold is exceeded, thus approximating larger perturbations, will 
be used to evaluate the designs. 
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Autopllot/Alrcraft S 1 mu latlon 


In order to verify the discrete control laws in a more complete 
fashion than Is possible with a linear analysis, a simulation of the 
aircraft/autopilot has been coded. It will enable us to determine the 
effect of the various autopilot saturations and the spoiler deadband 
on the autopilot performance and will provide an Independent verification 
of the designs. Furthermore , flight path time histories in the presence of 
wind gusts will be obtained to evaluate accuracy of the various designs. 

The simulation is a Fortran program of the aircraft controlled by 
the autopilot. The aircraft dynamic equations are the linearized small 
perturbation equations for the glide slope configuration written in 
standard first order form: 


x=Ax+Bu+B'w 

These are integrated by a fourth order Punga Kutta algorithm using an 
input initial state and time increment. The state vector expresses the 
vehicle’s angular velocity and orientation, as well as its rectilinear 
motion, and position in a ground fixed reference frame. A general auto- 
pilot subroutine determines the control vector on the basis of the current 
state after a user-specified number of repetitions through the integration 
loop (effectively determining sampling rate) . 

An optional digital simulation of atmospheric turbulenc is included 
and the effect of uniform gusts and gust distributions over the aircraft 
is modeled by B' in the equation above. 
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The gust vector w represents a disturbance to the system consisting 

of wind gusts (u , v w ) together with gust angular velocities (p , 

K 8 1 K g 

r g) • U g* v g , w g , and p^ are uncorrelated turbulence components 

while q and r are related to and ^ V B respectively. The six 
8 g 6x 6 X 

turbulence i mponents have a normal magnitude distribution with zero 
mean and a variance related to the t.ean wind, a function of altitude. 

The power spectral densities of the turbulence components are Dryden 
Spectra, which assume the longitudinal spatial autocorrelation function 
of gust veloc? ty to be of the form 

_ r 

R(r) > c^e ^ 

Here, o is the gust variance, and r is a general spatial coordinate. L, 
the characteristic scale length of the turbulence, is related in general 
to altitude. 

The continuous domain power spectral densities implied above (2,3) are 
digitally recreated by passing a Gaussian random number sequence through 
appropriate difference equations (4). 

The program is completely coded and in the final stages of checkout. 
Simulations have been made of an aircraft trajectory on a glide slope and 
found to give accurate results. Checkout of the system coupled to rvn 
autopilot and the r/ind model is now being checked out. 
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Roughness Evaluation 


The results of the discrete linear analysis of tne digitized s-plane 
designs show that the actuator plus short period aircraft dynamics give an 
accurate representation of the high frequency behavior of the full long- 
itudinal dynamics. Furthermore, all significant root movement due to 
sample rate is exhibitled by this simplified mode. The intent of the 
roughness evaluation is to determine the effect of slow sampling in 
terms of response roughness or Jerkiness as perceived by a pilot. These 
effects will be most noticeable in the high frequency modes of the air- 
craft. Therefore, the short period aircraft dynamics with the elevator 
lag has been selected for these studies. 

To more fully evaluate pilot response the high frequency approximation 
to the autoland and the control wheel steering modes will be evaluated. 

The autoland mode essentially places the pilot in a monitoring type of 
function and may produce different reactions compared to the control 
wheel steering where he will be required to execute various tasks. 

The holds to be considered are (Fig. 15) the standard zero order hold 
(ZOH) for a baseline, the first order hold, and the trir.ngular hold. The 
first order hold essentially requires two outputs representing the step 
in the control, £ u in Fig. 16, and the slope, k^. It shares the dis- 
advantage of the ZOH that step changes in the control occur. 
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Figure 15 

To eliminate this disadvantage, a triangular hold (TAH) is included 
which has the additional advantage of being particularly straightforward to 
mechanize, i.e. an integrator is placed between the Zou and the actuator. 

Figure 17 combines the different elements of the evaluation and 
indicates where the different lortions of the computations are to be carried 
out. The aircraft model is given by: 


q(s) 

6e(s) 


K(s+a) 

(s+b)(s +2CuO s+u?) 


where b = actuator lag (10 r/ sec) 

= short period dynamics (dependent on flight 

condition) 

B(s) 1_ 

q(s ) s 
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Figure 16 
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Analog Computers (EAI 231R's) 



Digital Computer (CDC 6600) 


Figure 17 
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The algorithms to be implemented for the holds are: 
ZOH: u(t) = u ;t A < t < t t +l 

FOH: u ( t ) = u 1 +k i (t-t t ); < t < t 


TAH: uCt) = ^ u^dt 

o 

The autopilot equations will basically be discretized versions 
of the pitch rate feedback through a washout (Fig. 6), proportional 
pitch attitude feedback, and the stick force input through a filter 
with associated control wheel steering logic. In addition, extra 
compensation for slow sampling will add one past value of measured 
states and one past value of control output for all holds, For the 
FOH, the control slope is required necessitating one additional 
difference equation: 


k 

i 


u - u 
i i-1 


T 


For the TAH, the additional integration of the plant will require 
further compensation in the autopilot. To eliminate effects of drift, 
we anticipate the additional feedback of the integrated control will be 
useful. 

Analytical evaluations of roughness have been proposed and 

are: 


RF. 


RF_ 


RF, 


RF. 
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Au 
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where the terms are defined in Fig. 16. These functions will be 


evaluated in the digital computer and correlated with the pilot 
evaluations of roughness. 

The exact nature of the evaluations with regards to initial 
conditions for the autoland, test maneuvers for the control wheel 
steering, and form of pilot evaluations, will be determined in 
collaboration with Langley personnel. 
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Financial Status 


Expenditures are slightly below projections and are 
indicated In Figure 18. Due to the higher Stanford expenditure 
rates through the summer months, we anticipate higher productivity 
through this period. 
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